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Abstract 

In this paper, asymptotic properties of the loss probability arc considered for an M/G/l/N 
queue with server vacations and exhaustive service discipline, denoted by an M/G/l/N -(V, 
E)-queue. Exact asymptotic rates of the loss probability are obtained for the cases in which 
the traffic intensity is smaller than, equal to and greater than one, respectively. When the 
vacation time is zero, the model considered degenerates to the standard M/G/l/N queue. 
For this standard queueing model, our analysis provides new or extended asymptotic results 
for the loss probability. In terms of the duality relationship between the M/G/l/N and 
GI/M/l/N queues, we also provide asymptotic properties for the standard GI/M/l/N model. 

Keywords: M/G/l/N queue, GI/M/l/N queue, serve vacations, invariant measure, loss 
probability, asymptotic behavior 
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1 Introduction 

For a queueing system with a finite system capacity N, the loss probability P\ oss {N), defined 
as the steady-state probability with which an arriving customer finds no free room in the queue 
(or buffer) and is forced to leave immediately, is an important performance measure in real 
applications, see e.g. [U [31 [TUl [16l [21] . Simple closed-form solutions for the loss probability are 
available only for a very limited number of special cases such as the M/M/c/K+c queue. For 
others, numerical computations and simulation methods are common tools for the analysis of 
the loss probability, but they cannot provide closed-form properties. Another important aspect 
is the asymptotic of the loss probability, which can lead to approximations, performance bounds 
among others. Also, approximations can be often obtained through asymptotic properties in 
the invariant measure for the corresponding infinite-buffer system, for which the analysis is 
usually less challenging than that for the finite-buffer one, see, e.g. [3j HJ [T2], [T7] . This paper, 
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investigating the asymptotics of the loss probability for an M/G/l/N queue with vacations, will 
go along this direction. 

Now we introduce this system based on the standard M/G/l/N queue. Recall that in the 
M/G/l/N queue, the arrival process is a Poisson process with rate A and the service times are 
i.i.d. random variables with a generic random variable S, which is independent of the arrival 
process, and has distribution function S(x) and probability density function s(x). The queue 
has a finite buffer of size N — 1 to store the incoming customs. Therefore, N is the maximal 
system capacity for holding customers including the one in service. To introduce the vacation 
model M/G/1/N-(V, E)-queue, we further assume in the M/G/l/N queue that whenever the 
system becomes empty, the server starts a vacation of random length V immediately. If, when 
the vacation ends, the system is still empty, the server takes another independent vacation of 
the same length as in the previous one until there is at least one customer waiting in the buffer 
upon return of the server from a vacation. From then, the server keeps serving customers until 
the system becomes empty (all customers have been served or exhausted). The vacation time 
V is assumed to be a non-negative random variable with the distribution function V(x). Then, 
p = XE[S] is the traffic intensity of the model. Throughout the paper, assume that E[S] < oo 
and E[V] < oo. Let Z + and C be the sets of the non- negative integers and the complex numbers, 
respectively. Denote by a,j and Uj the probability that j customers join the queue during a service 
time S and a vacation time V, respectively, i.e., for j 6 Z + , 



Various performance measures, such as the loss probability, the queue length, the waiting 
time, the busy period, the system throughput and so on, were investigated for the M/G/l/N- 
(V, E)-queue in the literature, e.g. [HI [UJ [HI [9]. The loss probability, which was studied in 
\13\ [TT[ [9], has been treated as a central issue among these performance measures. In these 
works, the loss probability for the finite M/G/1/N-(V, E)-queue was expressed in terms of the 
invariant probabilities for the embedded queue of the finite or infinite buffer queues, which are 
again unknown. As far as we know, for the model studied here no asymptotic analysis, as 
N — > oo, for the loss probability has been reported, which stimulated us to provide a detailed 
characterization of the asymptotic analysis of the loss probability. In the literature, asymptotic 
results for the M/G/l/N queue (without server vacations) can be found, for example, in [3]. 
In addition, asymptotic rates for the loss probability of the GI/M/l/N queue were obtained in 
(0 E] ) 5 by using the duality relationship between the M/G/l/N queue and the GI/M/l/N queue. 
However, asymptotic analysis of the loss probability seems not complete even for the M/G/l/N 
and GI/M/l/N queues and some important cases remain unexplored, which stimulated us to 
introduce new asymptotic methods to study the loss probability for the more general M/G/l/N- 
(V, E)-queue so that the methods or the results can be applied to derive new asymptotic behavior 
for those queues. 

In this paper we will provide a detailed characterization of asymptotic behavior of the loss 
probability for the M/G/l/N- (V,E)-queue. Some of the results are intuitively expected and 
others are not, and some of results are easy consequences from the existing literature and others 
deserve more detailed analysis. The characterization is obtained by investigating the dominant 
singular point in terms of singularity analysis, such as the Tauberian theorem, Tauberian-like 




and 
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theorem (Theorem A.l in Appendix) or a direct complex analysis. As special cases, for the 
standard M/G/l/N queue and the GI/M/l/N queue, we will provide new asymptotic results. 
Specifically, our main contributions include: (1) For p < 1 and under a light-tailed condition, 
two cases are studied in Theorem 13.11 for the first case, the decay is light but not exactly 
geometric, a new asymptotic behavior, and for the second case, decay is exactly geometric; (2) 
Also for p < 1, asymptotic results when either the service or the vacation variable is heavy-tailed 
are provided in Theorem 13.21 This case was not considered in [3]; (3) For p = 1, asymptotics 
are studied in Theorem 13.31 under the assumption that the second moment of the service time is 
finite and infinite, respectively. The infinite case was not considered in [3]; (4) An asymptotic 
result for p > 1 is reported in Theorem 13.41 which can be expected from the corresponding result 
in [3]; (5) When V = 0, the vacation model is degenerated to the standard M/G/l/N queue, for 
which asymptotic results are summarized in Theorem 14.11 We also provide asymptotic results 
for the GI/M/l/N queue in Theorem 14.21 by using the dual relationship with the M/G/l/N 
queue. These two theorems provide extensions of existing literature results. 

The rest of the paper is organized into: Section 2, which bridges the asymptotics of the loss 
probability of the finite queue and decay behavior of the invariant measure of the corresponding 
infinite queue; Section 3, which is devoted to asymptotic analysis of the M/G/l/N-(V,E)-queue 
and contains main results of the paper; Section 4, in which asymptotic results for the standard 
M/G/l/N queue the GI/M/l/N queue are provided; Section 5, which provides concluding re- 
marks; and an Appendix, which contains some properties in asymptotic analysis used in this 
paper for convenience. 



2 Connection with the M/G/l/oo-(V, E)-queue 



Let L(t) be the queue length of the M/G/l/N- (V,E)-queue at time t. Since the state space is 
finite, the system L(t) is always stable. It is known that L{t) itself is not a Markov process, but 
it can be analyzed through its embedded chain L{r n ) at the departure epochs r n , n £ Z+. The 
transition matrix of the embedded chain L(r n ) is given by 



P(N) 



b 




b 2 ■ 


bN-2 


1 


Of) 


a i 


a 2 ■ 


■ ClN-2 


1 





«o 


ai . 


■ O.N-3 


1 








a . 


■ 0,N-A 


1 








. 


a 





2^k=0 b k 
2^k=0 a k 
2_>fc=0 a k 
2^k=0 a k 



1 - ao / 



NxN 



where bj - ^ i=l 
probability invariant vector of P(N), i.e. tt(N)P(N) = tt(N), or 



Obviously, P(N) is positive recurrent. Let tt(N) be the 



3+1 



tt^N) = n (N)bj + K k (N)a j+l _ k , j = 0, . . . , N - 2, 



(2.1) 



k=l 



and E£o 1 ^(^) = l- 

Let 7T*{N) = limj^oo P[L(t) = j], j = 0, 1, . . . , N. By the Poisson Arrival See Time Average 
Property, we see that tt*(N) is also the probability that there are j customers in the system just 
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before an arrival, which implies that P\ oss (N) = tt* n (N). By Theorem 4.2 in [9], we have 

(1 - ^o)A- 1 



1 



(2.2) 



Therefore, the asymptotic analysis of the loss probability solely depends on the property of 
ttq(N). In the following we establish a relationship between the stationary probabilities for the 
finite-buffer queue and the invariant measure for the corresponding infinite-buffer queue for all 
values of the traffic intensity p. 

To connect the loss probability P\ oss (N) with the invariant measure of the corresponding 
infinite queue, we now introduce some related notations. For sequences {x n ,n G Z + } and 
{y n ,n € Z + } of non-negative real numbers, define 



S x {n) 



n-l 

i=o 



S x (n) 



oo 
j=n 



n > 1, 



and write x n ~ y n to imply limn._5.0o x n /y n = 1. Similarly, for two functions f(x) and g(x), 



f(x) ~ g{x) as x — )• xq— stands for lim 



X — >XQ- 



f{x)/g(x) = 1, and f(x) ~ g(x) for lim 2 



¥+00 



f(x)/g(x) 



1. For any distribution function F(x) on the nonnegative real line, let f(x) be its probability 
density function when it exists, F(x) = 1 — F(x) be its complementary distribution function, 



F*(t) = 

be its Laplace transform, and 



-tx 



: dF(x) 



-tx 



f(x)dx, t G C 



F e (x) 



E[X] 



F{u)du 



be its equilibrium distribution whenever E[X] is finite. The distribution function F(x) is said 
to be regularly varying if F(x) = x~ a L(x), where L{x) is a slowly varying function that satisfies 
lim. r _, +00 = 1 for any t > 0. 



The corresponding infinite-buffer vacation model is referred to as the M/G/l/oo-(V, E)- 
queue, and its embedded chain has the following transition matrix: 



P 



where the coefficients 6j and are the same as those in P(N). Regardless of stability or non- 
stability of the infinite-buffer queue, P has a unique invariant measure, up to a multiplicative 
constant, due to its special structure. Let ir be the invariant measure of P, i.e. irP = ir, or 
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(2.3) 



Define the generating functions 
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for bk, ak and tt^, respectively. It is routine to calculate 

A W = S '(A-Az), BM= S-(A-A--)|V(A-A,)-, 1 

2(1 - ^o) 

From (12. 3p . we have 



n(z)[z-A(z)]=7T [zB(z)-A(z)}, 
whenever H(z) is finite, that is, 

U(z)[z - S%\ - \z)\ = MV*(*-^)-l]S*(\-\z)^ 

1-Uq 



Observing that all equations of (J2JJ) are identical to the corresponding ones in (|2.3p . we have 

7r *( JV ) = ^7M' * = 0,1,---,JV-1. (2-5) 

Remark 2.1. iVote i/iai P is positive recurrent, or ir is a probability measure, if and only if 
p<l. 

The following lemma will later be frequently used to analyze the asymptotic behavior of the 
loss probability in Section 3. 

Lemma 2.1. For the case of p < 1, we have 

lim - 1 P] n ^( N ) = 1 - P; ( 2 -6) 
N^ooS*{N) loss 

and for the case of p > 1 , we have 

E[V]tt 



lim SJN) 

JV-»oo 



*w*o - ( 1 



pE[S](l-u ) 



(2.7) 



Proof. Inserting (2.4) into (2.2) yields 



p , m _ E\V]no + (E[S] - A^Kl - v )S n (N) 
lossl 1 E[V]^ + E[S]{l-v,)S^N) ' { ^ 

When p < 1 , we have X)£o n i = ^ anc ^ = aH [ V ] ^"^ " ^ se Q2.8p to derive 

p , N) = d-P)SAN) 

from which (|2.6p follows easily. 



When p > 1, we have from (|2.8p 

lim P loss (iV) = 1 - -, 
since S 7r (N) -> oo as iV ->■ oo. Write ([ZgJ) as 

£[y]vr 



^loss(^) - (l - 



p^[V]7ro + p5 ff (JV)^[5](l-m)' 
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Then, (|2.7p follows immediately. □ 

From Lemma 12.1^ we see that the asymptotic behavior of the loss probability P[ oss (N) as 
TV" — > oo relies on the decay property in the invariant measure tt. For a stable system, we study 
the decay rate of the tail probability S n (N), and for a non-stable system, we analyze the decay 
of l/S n (N). 

For the Laplace transform F*(X - Xz) (F*(X - Xz) = F*(t)\ t =\-\ z ), define Rp* > 1 to be 
the leftmost singular point of the functions F*(X — Xz), 

F*'(X - Xz) = dF * (X - Xu) \ u=z = r Xte~^f{t)dt, 
du J 

and 

F*"(X - Xz) = d2F *fj Xu) \u= z = J o °°(Xtfe-^f(t)dt. 

Observe that S* (A — A • 1) = 1, 5*'(A — A ■ 1) = p and S*(X — Xz) is a strictly increasing and 
convex function about z when z > 0. This observation yields the following lemma, which will 
later be used to investigate the decay of S W (N) and l/S 7r (N). 

Lemma 2.2. (i) If p < 1 and 1 < Rs* < S*(X — XRs*), then there exists at most one solution 
z\ to the equation z = S*(X — Xz), 1 < z < Rs*- 

(ii) If p > 1, then there exists exactly one solution Z2 to the equation z = S*(X — Xz), 
< z < 1. 



3 Asymptotic behavior for the M/G/1/N-(V, E)-queue 

We perform the asymptotic analysis according to three cases: (i) p < 1; (ii) p = 1; and (iii) 
p>\. 



3.1 The case when p < 1 

For this case, we investigate the asymptotic behavior for two subcases according to min{i?s* , Ry* } 
> 1 or m.in{Rs* , Ry* } = 1- The former subcase corresponds to light-tailed characterizations, 
for which the asymptotic rates are either exactly geometric or not exactly geometric. The latter 
subcase corresponds to that one of the sequences {b{\ and {aj} is heavy-tailed, for which a 
regular varying condition is imposed on the service time and vacation time distributions. 

Theorem 3.1. Suppose that p < 1 and r := min{i?5* , Ry* } > 1. 



(i) Suppose that r = Ry* and S*(X — Xr) < r. If V*(X — Xz) can be analytically extended 
to a A(r)-domain (see Definition A.l in Appendix) and for some 8 > 0, V*(X — Xz) ~ 

( r - z ) s ' as z ~ ^ r ' z e ^( r )> then 

hmiV^P, (N)= c(l- P ) 2 S*(X-Xr) 

urn^iy Hoss^l r e ~ 1 T(8)XE[V](r — l)[r — S*(X — Xr)} ' { } 
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(ii) Suppose that r = R s * and r < S*(X — Xr), then there exists exactly one solution z\ to the 
equation z = S* (A — Xr), 1 < z < r, such that 



lim z N Pi (N) - *i(l-P) 2 [l-^-A.i)] 
lun^z, P l0SS {N)- XE[v]{zi _ 1)[1 _ s *, {x _ Xzi) y 



Proof, (i) Due to the assumption, we know that V*(X — Xz) is analytic in a A(r)-domain, which 
implies that there exists some r±, T\ > r and some a, < a < \ such that V*(X — Xz) is analytic 
in the domain A(r, ct,r\). Thus we can find an r 2 such that r < r 2 < r\ and V*(X — Xz) is 
analytic in the domain A(r, a,r 2 ). Hence by the expression (|2.4p . we know that H(z) can be 
analytically extended to this domain A(r, a,^), i.e. H(z) is A(r)-analytic. 

Since V*(X — Xz) ~ , _f , fl as z —> r in the complex domain A(r, a, r 2 ), we have 
limIL(z)(r - z)° - ^"^(A-Ar) 



XE[V][r - S*(X- Xr)} 



that is, 



E n/ z N „ c(l - p)S*(X - Xr) 1 . , . 

^ ( r } ~ ASM[r-g«(A-Ar)]r»(l-*)* ' aS ^^^(r,a,r 2 ). 



n=0 



It follows essentially from a result in [8] (see Theorem lA.il in Appendix) that 

„ „ c(l-p)S-(A-Ar) 

n r fi r(fl)A£?[Vl[r-5*(A-Ar)] 
from which and the Stolz-Cesaro theorem, we have 

S (n) ~ c(l-p)g*(A-Ar) e _ x 

Then, assertion (|3.ip follows from Lemma 12.11 immediately. 

(ii) It follows from Lemma 2.2 that there exists exactly one solution z\ to the equation 
z = S*(X - Xz), 1 < z < r. Since 

. zi(l-p)[l-V*(X-Xz 1 )} 

lim U(z)(z 1 -z)= v F ;/ L 

VA ; AE[F][1-5*'(A-A«i)] 

zi is a simple pole for II(z). For any z 6 C such that |z| = zi and z / zi, we have 

/>oo 

|5*(A - As) | < / e- Ax |e A2a; |^) < S*(X - X\z\) = z u 
Jo 

which implies that there is only one solution z = z\ to the equation z = S* (X — Xz) on the circle 
\z\ = z\. Hence z = z\ is the only singular point of H(z) on the circle \z\ = Z\. By Theorem 5.2.1 
in [20], a standard result on asymptotics of complex functions, we have 

(l-p)[l-V*(X-Xz 1 )] 
n ~ A£7[y][l-5*'(A-M)] 1 ' 
from which and the Stolz-Cesaro theorem, it follows that 

Tj / x zi(l-p)[l-F*(A-Azi)] _ n 

a^v^i-i^i-s^a-a^)] 1 ■ 

The assertion (13. 2D follows from Lemma 12 . 1 1 immediately. □ 
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Remark 3.1. (i) In the above theorem, the assumption S*(X — Xr) < r in (i), which is equivalent 
to that there is no solution to the equation z = S* (A — Xz) for 1 < z < r, ensures that the 
dominant singular point of 11(2;) is determined by V*(X — Xz), while the assumption in (ii) 
ensures that the dominant singular point ofH(z) is determined by the solution to the equation z = 
S*(X — Xz) for 1 < z < r. (ii) In the first case, the exact asymptotic behavior of Pi 0SS (N) reveals 
a non-geometric phenomenon when 6^1, which is faster than the geometric rate r~ N when 
9 < 1 and slower than the geometric rate r~ N when 6 > 1. This is a different characterization 
from the well-known exact geometric decay in the related literature. 

Example 3.1. We now give an example to illustrate that how to verify the analytic condition 
in (i) of Theorem 3.1. Assume that (i): X = 1, V has probability density function v(x) = 
~ k\ X e ~ ax -> a > 0, x > 0, and S has probability density function s(x) = pe~ px ,p > 0, x > 0. 
Then it is easy to calculate 

E\V] = —, E[S} = ~, P = XE[S] = l/p, 
a p 

V*(X - Xz) = V*(l -z)= e-^ x v(x)dx = — Re(z) < 1 + a, 

Jo (1 + a - z) k+L 

S*(X-Xz) = S*(l-z)= [°° e-( 1 - z >s(x)dx= ? , Re(z) < 1 + p. 

Jo 1+p-z 

So, we have Ry* = 1 + a, Rs* = l + p. Assume that (ii): p > 1 + a, then we have r = Ry* and 

S*(X — Xr) < r. Obviously, the function ^ s analytic m C \ {1 + &}, i.e. the whole 

complex plane except for the point z = 1 + a. Hence, the Laplace transform V*(X — Xz) can be 

analytically extended to C\{1 + a} , and the A(r)-analytic condition holds automatically. Thus, 

under the assumptions (i) and (ii), every condition in (i) of Theorem 3.1 is satisfied, so we have 

lim N- k (l + a) N P ]ns JN) = - ~ ) /p)2 -. 

JV^oo v ; lossK ' (l + a) k (k + l)\{p-l-a) 



and 



□ 



Theorem 3.2. Assume that p < 1 and min{i?s* , Ry* } = 1. 
(i) IfV{x) ~ cS(x) ~ x~ a L(x) for some c > and a > 1, then 



(ii) IfV{x) = o(S(x)) and S(x) ~ x a L{x) for some a > 1, then 



N a-1 

lim — —p, (N) = . (3.4) 

tv^oo L(f ) loss a - 1 y J 



(iii) If S(x) = o(V(x)) and V(x) ~ x a L(x) for some a > 1, then 



lim (N) = ± '{ r - . (3.5) 

N^oo L(f ) lossy 1 (a - 1)E[V] K ' 
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Proof, (i) It was shown in Proposition 1.5.10 in [5] that if F{x) ~ x a L(x) for some a > 1, then 
the equilibrium distribution 

— . , x-^-V L(x) 

Using L'Hospital rule obtains 

lim ZM - hm m(E[V]-I X V(t)dt) _ E[S] 

5 e (x) sry](£?[5] - /* 5(t)«te) «-^oo 5( x ) ' v • ; 

whenever the limit Hm x _>. C i exists. Using (j3.6[) and (j3.7[) yields 

U e (x) ~ S e (x) ~ 



E[V] ev ' E[V](a-l)' 
Then by Proposition 6.4 in [2], we have 

S (k) ~ f— + g ^ >L±k-<-~-*>L ( k - 
bn[k) \E[V] + c(l-p)E[S]J a- I L {\ 



from which and Lemma \2.1\ (|3.3p follows, 
(ii) It follows from (JHSJ) and ((221) that 



Se(x) ~ ^fgTT^ _ . V e (a:) = o(S e (x)). 



From Proposition 6.1 in [2], we derive 



from which and Lemma 12. 11 (13. 4p follows, 
(iii) It follows from and ((277]) that 

From Proposition 6.2 in [2], we derive 

from which and Lemma 12,1} (|3.5[) follows. □ 

Remark 3.2. (%) Under the assumption that S(x) ~ x~ a L(x) (V(x) ~ x~ a L(x)), it is natural 
to restrict that a > 1 due to E[S] < oo fE[U] < oo), the basic assumptions of this model, (ii) 
Heavy-tailed asymptotic results were obtained for M/G/l-type Markov chains (e.g. fl5\. [TP ) /). 
in terms of asymptotics of a(n) and b(n), which can not trivially lead to asymptotic results in 
terms of queueing primitives or service/ arrival parameters, for example, results presented in this 
paper. 
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3.2 The case when p = 1 



For this case, if the second moment of the service time is finite, we provide the asymptotic decay 
for the loss probability. When the second moment is infinite, we consider a family of service 
distributions to understand the asymptotic behavior. 

Theorem 3.3. Assume that p = 1. 

(i) IfE[S 2 ] < oo, then 

^ NP loss( N ) = (3-8) 

(ii) If E[S 2 ] = oo, let s(x) ~ cx~ 8 for 2 < 9 < 3, where c is a constant. Then, for 2 < 9 < 3, 

iToc^ P lossW- (i_0)(2- 0)(3-6l) ' (3 ' 9) 

and /or 9 = 3, 

N c\ 2 
lim — -P /ns „(AT) = i— . (3.10) 

Proof, (i) Since 

z-S*(X-Xz) , 1-5*'(A-A^) -S*"(X — A • 1) A 2 MS 2 ] 

lim ; -7s = Inn = = , 

(l-^) 2 -2(1 -z) 2 2 

we have 

lim nw(1 _ _ . lim "*<* -*>-i- 2 *»^1 



wv ; (1 - ZA>)A 2 £[S 2 ] *-n- 1 - z (1 - z, )A 2 £[,S 2 ] ' 

It follows from a Tauberian theorem (Theorem IA.3I in the Appendix) that 

2^XE[V]n 
SAn) ~ (1-, )A 2 ^ 2 ]- 

From Lemma |2.H we have 

which is (J3IHD by noting that p = XE[S] = 1. 

(ii) Let C/(x) = f*(Xt) 2 s(t)dt. There are two cases. 
(1) If 2 < 9 < 3, then we have 

U(x) cX 2 



lim 



x->oo x 3 s 3 — 9 
Obviously, 

roc roc 

S*"(X-Xz) = / (Xt) 2 e- ( - x - Xz ^ t s(t)dt= / e-^ Xz ^dU(t). 
Jo Jo 
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Using a Tauberian theorem (Theorem IA.2I in the Appendix) yields 

By L'Hospital rule, we have 

z-S*(X-Xz) , l-S*'(X-Xz) -cA 9 " 1 r(4 

lim — — 7- — = hm 



(1 ^l-_(l_0)(l_^)-(2-0) (l-0)(2-0)(3 
from which it follows that 



(2 _ e) [^( A -Az)-l](l-zr( 2 - e ) 

(l-i/ ) z "i- -cA«-ir(4-e) 
v y (i-e)(2-e)(3-e) 

tt o (1-0)(2-0)(3-0)A£;[V] 



lim n z (1-2 )-V-°> = - lim 

y ' (1-0) (2-0) (3-0) ^ Z J 



c(l - i/o)A 9 " 1 r(4 - 0) 
Now, by using another Tauberian theorem (Theorem IA.3|) . we have 

ir (l-9)(2-6)(3-6)XE[V] n e ~ 2 



S n (n) 



c(l - i/ )A e ~ 1 r(4 - 6») r(0-l)' 
from which and Lemma \2.1\ (|3.9p follows immediately. 
(2) If 9 = 3, then 

hm — ^ = cA . 
Using a Tauberian theorem (Theorem IA.2[) yields 

S*"(X - Xz) ~ cA 2 In f A _^^ = -cA 2 ln(A - Az). 
Then, by L'Hospital rule, we obtain 



and 



lim n*-A»)-i = 

Z-l 



z-S*(A-Az) , 1-S*'(A-Az) 

hm -77- — 7- — r = hm 



(l-z) 2 ln(A - Xz) (1 - z) [2 ln(A - Az) - (1 - z)] 

lim —S*"(X — Az) 

~~ 21n(A - Az) + 2(1 - z) + 2 
cA 2 
2 



It follows from (I3TTTD and (I3TT21 that 



(3.12) 



27T £[U] 

c(l - f )A' 
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Hence, we have 



TTf ^ 1 2tt E[V] 1 

{Z) ~{l-z) Ht&s) ~ " -o)A (1 - z) ln(^) ' &S * 

Using another Tauberian theorem (Theorem IA.3|) leads to 

2tt q E[V] n 



SJn) 



c(l — uo)X Inn 

from which and Lemma |2. 11 (|3.9p follows immediately. □ 

Remark 3.3. (a) In the assumption of s(x) ~ cx^ 61 in case fiij, if is reasonable to restrict 9 
in the domain (2,3]. Otherwise, if 9 = 2, then E[S] = oo, and if 9 > 3, then E[S 2 ] < oo, both 
contradict other assumptions, (b) The condition that E[S 2 ] < oo is often imposed in studies of 
this type of problems, for example in J3}/, but as far as we know, Case (ii) is a new one. (c) 
In the case p = 1, the asymptotic behavior of the loss probability is independent of the vacation 
time V . 



3.3 The case when p > 1 

This is the simplest case, in which we no longer need to impose any further property assumptions 
on the service time or on the vacation time for the following asymptotic property for the loss 
probability. 

Theorem 3.4. Assume that p > 1. Then we have 

(l-z 2 )[l-S*'(X-Xz 2 )]E[V} 



lim Zr, N 

N->-oo 



Ploss( N ) - (l - £) 



z 2P E[S][l-V*(X- Xz 2 )} 
where z 2 G (0, 1) is the unique solution to the equation z = S* (A — Xz). 



(3.13) 



Proof. Prom Lemma 2.2, we know that there exists a unique solution z 2 E (0, 1) to the equation 
z = S* (X — Xz). Obviously, S*'(X — Xz 2 ) < 1. Since 

r tt/ w s z 2 vr [l-y*(A-Az 2 )]] 
hm IL(z){z 2 -z ) = - — —j- — T , 

z^Z2- (1 - U )[l — S* (X - XZ2)\ 

z 2 is a simple pole for n(^r). Similar to the proof of Theorem 13. 1\ we know that z = z 2 is the 
unique singular point of H(z) on the circle \z\ = z 2 . By Theorem 5.2.1 in |20| . a standard result 
on asymptotics of complex functions, we have 

7ro[l-V*(A-A*2)] _ n 

vr„, ~ 



(l-u )[l-S*'(X-Xz 2 )} 2 ' 
from which and the Stolz-Cesaro theorem it follows that 

q ^ z 2 t: [1-V*{X-Xz 2 )} 

b ^ n >~ (l- Z2 )(l-u )[l-S*'(X-Xz 2 )f 2 ■ 

The assertion fj3. 13H follows mow from Lemma 12.11 immediately. □ 
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4 Decay for standard M/G/l/N and GI/M/l/N queues 



The M/G/l/N queue with vacations becomes the (standard) M/G/l/N queue when the vacation 
time V = 0. The embedded transition matrix of the M/G/l/N queue is obtained by replacing 
bk by afc, k = 0, • • • , TV — 2, in the first row of the embedded transition matrix for the M/G/l/N 
queue with vacations. Denote by P\ oss (K) the loss probability for the M/G/l/N queue. It is 
known that 

where n(N) is the invariant vector of the embedded Markov chain for the M/G/l/N queue. 
The following result is obtained by the same asymptotic analysis for P\ oss (N) in Section 2 with 
V = 0. 

Theorem 4.1. For the M/G/l/N queue, we have the following asymptotic results on the loss 
probability (iV). 

(i) Let p < 1. 

(1) If Rs* = 1 and S(x) ~ x~ a L(x) for some a > 1, iaen 



^T7NT P loss( N ) 



TV^ooL(f) ' a-1 



(2) If Rs* > 1 and Rs* < S*(X — XRs*), then there exists only one solution o~\ to the 
equation z = S* (A — Xz), 1 < z < Rs* such that 



lim o-?P, ns JN) = — ^1- 4 ir 
tv->oo 1 tossV ; [S* (X - Ao-i) - 1] 



(ii) £e£ p = 1. 

(%) //P[S 2 ] < oo ; then 



^ Np loss( N ) 



X 2 E\S 2 } 



N^oo ' 2 

(2) If E[S 2 ] = oo and s(x) ~ cx~ e , iaen for2<9 <3, 

hmiv^A (JV) = ^ 1 r(g-i)r(4-g) 
amV tesl j (i-e)(2-e)(3-e) : 



and /or = 3, 



(iii) Lei p > 1. f/ien 



JV A cA 2 
iV-K» In 



lim cr 9 

JV-)-oo ^ 



V 



PlossW - (l - J) 



1 - S*'(A - Aa 2 ) 



02P 2 



where o~2 € (0, 1) is i/te unique solution to the equation z = S*(X — Xz). 
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Remark 4.1. (a) The above results for p > 1 and for p = 1 urai/i the assumption that E[S 2 ] < oo 
coincide with those in (JEj). (b) When p < 1, we show that the asymptotic rate can be exactly 
geometric or polynomial (heavy-tailed). The later case was neglected in (c) The asymptotic 
result for the case that p = 1 with E[S 2 \ = oo is considered new. 

We now consider the (standard) GI/M/l/N queue. Denote by A the generic interarrival 
time random variable, and by A(x) and a{x) its distribution function and probability density 
function, respectively. Let p be the service rate and let N be the maximal capacity of the 
system, or the number of customers in the system including the one in service. The traffic 
intensity of the GI/M/l/N queues is p = ^[A ] • ^ loss^) we denote the loss probability of 
the GI/M/l/N queue. Using the dual concept (see, e.g. (3J HE]), we have 



Pl oss (N) = n N (N) = 7r (N + l), 



(4.2) 



where tt(N) is the invariant distribution of the embedded GI/M/l/N queue, and ft(N + 1) is 
the invariant distribution of the M/G/l/N + 1 queue in which the arrival process is Poissonian 
with parameter p, the distribution of the service time A(x) and the traffic intensity p = 
From (12. 5ft . (14. 2j) and the proof of Theorems l3.2fj3T3| we obtain the following assertion. 

Theorem 4.2. Consider the GI/M/l/N queue. 



(i) Let p < 1. We then have 



lim ViPlossW = 1 " A (P ~ UK) 

iv— >oo 



where rj\ E (0, 1) is the unique solution to the equation z = A*(p — pz). 

(ii) Let p = 1. 

(1) If E[A 2 ] < oo, then 



lim NP loss {N) 



p 2 E[A 2 



(2) If E[A 2 ] = oo and a(x) ~ cx~ e , then for 2 < 9 < 3, 

cn 9 - x Y{e - i)r(4 



^N»- 2 P los JN) 



JV->oo 



(l-0)(2-0)(3-0) 



and for 



lim JLpj (N) = 

TV^oolniV Lossy ' 2 



(4.3) 



(iii) Let p > 1 and Ze£ i?^* 6e i/ie leftmost singular point of the function A(p — pz). 
(1) If Ra* = 1 and A(x) ~ x~ a L(x) for some a > 1 as x — >• oo, £/ien 



2V 



a-l 

lim 

AT^oo L(f ) 



A a 
a — 1 



(2) If Ra* > 1 -Ra* < 5'*(A — Af?^*), then there exists only one solution 772 to the 
equation z = A*(p — pz), z € (1,Ra*) such that 



lim 7? 9 



N 



P loss( N ) - -p 



A*(p - pn 2 ) - 1' 
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Remark 4.2. (a) The results for p < 1 and for p = 1 with E[A 2 ] < oo coincide with those 
in J3J EJ/. (^J When p > 1, we showed that the asymptotic rate can be either geometric or 
polynomial. The latter case was not considered in (JBj). (c) The asymptotic result on the case 
that p = l with E[A 2 ] = oo is considered new. 



5 Conclusions 



The main contribution of this paper is a detailed characterization of the asymptotic rates of the 
loss probability for the M/G/l/N-(V,E)-queue. New asymptotic behavior was found even for 
simpler M/G/l/N queue and the GI/M/l/N queue in terms of introduction of new asymptotic 
analysis methods. 

Although we only focused on asymptotics of the loss probability, the arguments in this paper 
can be extended to derive asymptotic results of the other performance measures, such as the 
queue length distribution. Let L(t) be the queue length process of the M/G/l/oo-(V, E)-queue. 
Suppose that L(t) is stable, i.e. p < 1. It is well known that L(t) has the invariant distribution 
7r given by (2.3), i.e. L(t) has the same invariant distribution as that of its embedded queue. 
Define 

N oo 

h*(N)-*\\=Y,\KH N )-«i\+ E 

3=0 j=N+l 

where tt^(N) is given by (2.2) and tt*(N) is given by (see formula (13) in [9]) 
n*(N)= ^(I-HQA- 1 <i<iV-l 

Since TVj(N) = s n ^ N ^ , < j < N — 1, and ttq = xe[V]' ^ ' ^ ^ s no ^ nar d to find that the 
asymptotic rate of ||7r*(iV) — 7r|| as N — > oo is completely determined by the decay rate of S n (N) 
as N — > oo. Hence the asymptotic rate of ||7r*(iV) — tt|| can be derived by using similar arguments 
to that in Section 3.1. 
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A Appendix 



The following Definition A.l and Theorem A.l are taken from Definition VI. 1 and Corollary 
VI. 1 in [8], respectively, with a slight difference. Theorem A.l is called Tauberian-like theorem 
in 



Definition A.l. Let r > 1. For given constants (ft, f with r > r and < /3 < S, define A(r, /3, r) 
to be the following open domain 

A(r, /?, f) = {z : \z\ < f, z ^ r, \Arg(z — r)\ > f3}. 

A domain is called a A(r)-domain if it is a A(r, /3,f) for some f and (3. A function is said to 
be A(r)- analytic if it is analytic in some A(r)-domain. 

Theorem A.l. Assume that f(z) is A(l)- analytic and 

/0) ~ 7: ° 7g , as z->l, zeA(l) 

for some constant c and 9. Let f n be the nth Taylor coefficient of f(z). (i) If 6 ^ {0, —1, —2, • • • }, 
then 

f c 0-1 

fn- W) n ; 

(ii) If Be {0, -1,-2, • • • }, then 

fn = o(n 6 - 1 ). 

The following theorems, referred to as Tauberian theorems, are taken from Section 5 of 
Chapter 13 in [7]. 

Theorem A. 2. Let U be a measure on [0, oo] with Laplace transform uj(t) = e~ Xx U{dx}. 
If L is a slowly varying function and < p < oo, then the properties 

w(r) ~ t -pl(-\ , r^O, 



and 

U(t) ~ ^-— tPL(t), t^oo, 

are equivalent. 

Theorem A. 3. Let q n > and assume that 

oo 

Q( S ) = Y,QnS n 
n=0 

converges for < s < 1. If L is a slowly varying function and < p < oo, then the properties 
and 



+ qi-\ h q n -i ~ w 1 , ,, n p L(n), n ^ oo, 

r(p + i) 



are equivalent. 
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